Collective Modes and the Superconducting State Spectral Function of Bi2212 
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Photoemission spectra of the high temperature superconductor Bi2212 near (n, 0) show a dramatic change 
when cooling below T c : the broad peak in the normal state turns into a sharp low energy peak followed by 
a higher binding energy hump. Recent experiments find that this low energy peak persists over a significant 
range in momentum space. We show in this paper that these data are well described by a simple model of 
electrons interacting with a collective mode which appears only below T c . 

PACS numbers: 71.25.Hc, 74.25.Jb, 74.72.Hs, 79.60. Bm 
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Angle resolved photoemission spectroscopy (ARPES) 
has become one of the key tools used to elucidate the 
physics of high temperature superconductors. It has pro- 
duced a number of important observations concerning the 
nature of the normal and superconducting phases. Ex- 
amples are the existence of a large Fermi surface, and 
an anisotropic enpjrgy gap in the superconducting and 
pseudogap phaseslil. The most interesting aspect of the 
ARPES data, though, is the unusual nature of the spec- 
tral lineshape and how this lineshape changes as a func- 
tion of doping, momentum, and temperature. Perhaps 
the most profound example in this regard is the temper- 
ature dependence of the lineshape near the (tt, 0) point 
in Bi2212. A very broad normal state spectrum evolves 
quite rapidly below T c into a resolution limited quasi- 
particle peak, followed at higher binding energies by a 
dip then a hump, after which-the spectrum is equivalent 
to that in the normal stateJjQ Similar effects have been 
seen in tunneling spectra, where it has been found that 
all of these spectral features (peak, dip, hump) scale with 
the superconducting gapO. This implies that the electron 
self-energy has a dramatic change below T c . 

In a recent paper, our group has shown that the low en- 
ergy peak persists over a surprisingly large range in mo- 
mentum space along the (tt, 0) — (0, 0) and (tt, 0) — (tt, tt) 
directionso. As argued in that paper, this result can be 
connected to the change in lineshape with temperature 
noted above. The idea is that the dip in the spectrum at 
(tt, 0) implies that the imaginary part of the self-energy, 
IrnE, has a step-like drop from a large value at binding 
energies larger than the dip to a small value for smaller 
energies. This step behavior has recently been verified 
by us by a direct extraction of S from the dataQ. By 
Kramers-Kronig transformation, then, ReT, will have a 
strong peak at the dip energy. The consequence of this 
is that there will always be a low energy quasiparticle 
pole trapped on the lower binding energy side of the dip 
energy, even when the normal state binding energy is 
quite large. It is this effect which we believe leads to the 
persistent peak. 

This raises the question of what kind of microscopic 
picture can lead to such behavior. As discussed in our 



paperi, a step edge in JmS is equivalent to the prob- 
lem of an electron interacting with a sharp (dispersion- 
less) mode. This model has been worked out in detail in 
the classic literature of strong-coupling^uperconductors, 
where the mode is an Einstein phonond. In the current 
case, though, the effect of the mode only appears below 
T c , and therefore implies a collective mode of electronic 
origin. Still, the mathematics is largely equivalent. What 
we show in this paper is that this simple model gives a 
good quantitative fit to the data. 
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FIG. 1. a 2 F for three models, MFL (dashed line), gapped 
MFL (dotted line), and gapped MFL plus mode (dotted line 
plus 5 function). Inset: Feynman diagram for the lowest order 
contribution to E from electron-electron scattering. 

We begin by discussing self-energy effects in supercon- 
ductors. For now, we ignore the complication of mo- 
mentum dependence. The lowest order contribution to 
electron-electron scattering is represented by the Feyn- 
man diagram shown in the inset of Fig. 1. In the su- 
perconducting state, each internal line will be gapped 
by A. This implies 4hat the scattering will be sup- 
pressed for |cj| < 3A0. This explains the presence of 
a sharp resolution-limited quasiparticle peak at low tem- 
peratures. What is not so obvious is whether this in 
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addition explains the strong spectral dip. —Explicit cal- 
culations show only a weak dip-like featurecll To under- 
stand this in detail, we equate the bubble plus interaction 
lines (Fig. 1 inset V,to an u a 2 F" as in standard strong- 
coupling literatureB. In a marginal Fermi liquid (MFL) 
at T=0, a 2 F(Q) is simply a constant in f2. The effect of 
the gap is to force a 2 F to zero for < 2A. The question 
then arises where the gapped weight goes. It could be 
distributed to higher energies, but in light of the above 
discussion, we might expect it to appear as a collective 
mode inside of the 2A gap. For instance, this indeed oc- 
curs in FLEX calculations where the bubble represents 
spin fluctuationsEj, in which case a sharp mode will ap- 
pear if the condition 1 — Uxoi^i^) = is satisfied for 
fl < 2A. These three cases (MFL, gapped MFL, gapped 
MFL plus mode) are illustrated in Fig. 1. 
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FIG. 2. (a) ImE for MFL (solid line), gapped MFL 
(dotted line), gapped MFL plus mode (dashed line), and 
simple d-wave model (dashed-dotted line). Parameters are 
a=l, w c =200meV, A=30meV (0 for MFL), fi =2A, and 
ro=30meV. (b) Spectral functions (times a Fermi function 
with T—14K convolved with a resolution gaussian of a— 7.5 
meV) for these four cases (e=-34meV). 



E is easy to obtain analyticallyu if we ignore the com- 
plication of the superconducting density of states from 
the k — q line of Fig. I and just replace this by a step 
function at A. The resulting TmE for the gapped MFL 
and gapped MFL plus mode models are shown in Fig. 2a 
in comparison to the normal state MFL. Note that struc- 
ture in a 2 F at SI appears in E at \uj\ — il + A due to the 
gap in the k q line. Moreover, the MFL plus mode is 
simply the normal state MFL cut-off at 3A (this is ob- 
tained under the assumption that all the gapped weight 
in a 2 F shows up in the mode). In contrast, the gapped 
MFL decays linearly to zero at 3A. This simple picture is 
not changed that much if one actually solves the strong- 
coupling equations for £ and At£l. (For an s-wave gap, 
the linear behavior of TroE is replaced by a square root 
behavior). 

The spectral function is given bj 



with (a complex) Z(u>) = 1 — Y,(lu)/lu. These are shown 
in Fig. 2b and were convolved with a gaussian of er=7.5 
meV, typical of high resolution ARPES, with a constant 
ImE (T ) added for |w| > A to reduce the size of the 
quasiparticle peak. We note that there is no dip as such 
for the gapped MFL model, whereas the addition of the 
mode causes a significant dip. The latter behavior is 
consistent with experiment. Moreover, the mode model 
has the additional advantage that JmE recovers back to 
the normal state value by 3A, which is also in agreement 
with experiment in that the normaLand superconducting 
state spectra agree beyond 90 meVQ. 

We contrast this behavior with that expected for a sim- 
ple d-wave model. To a first approximation, this can be 
obtained by replacing the step drop in JmE in the MFL 
plus mode model with (\oj\ - A) 3 for \u\ < 3AE3. This 
is shown in Fig. 2a as well, with the resulting spectrum 
in Fig. 2b. Only a weak dip appears. Moreover, we have 
analyzed models with the exponent 3 replaced by some 
n and have found that n must be large to obtain a dip 
as strong as seen in experiment. Therefore, the upshot is 
that at the least, something similar to a step is required 
in JmE to be consistent with experiment. 

In principle, we could take the above MFL plus mode 
model and fit experiment with it. In this paper, though, 
we consider a simpler model. There are several rea- 
sons for this. First, the MFL model has a number of 
adjustable parameters associated with it. There is the 
coupling constant (a), the cut-off frequency (w c ), and 
the mode energy (which is not in general 2A). More- 
over, the spectrum for k points near the (tt, 0) point does 
not appear to be MFL-like in nature. We have found 
that the normal state Bi2212 spectrum is fit very well 
by a Lorentzian plus a constant in an energy range less 
than 0.5eV. This is also true for Bi2201 spectrum where 
the normal state can be accessed to much lower tem- 
peratures. The constant term represents the so-called 
"background" contribution, and is essentially equivalent 
to spectrum for k > kp, where it is also seen that the 
background gets gapped by A in the superconducting 
state. There are several possibilities for what the back- 
ground could be due to, and in fact could be a combina- 
tion of all of these: (I) incoherent part of A, (2) inelastic 
secondaries, (3) emission from the BiO layers, (4) diffrac- 
tion of the photoelectrons by the surface BiO layer, etc. 
Since this has little to do with the peak/dip/hump struc- 
ture, we choose to subtract this off, but note the caveat 
that this is an incomplete description if part of the back- 
ground is intrinsic. Finally, the Lorentzian simplification 
allows us to directly obtain the dispersion et fxom tight 
binding fits to the normal state peak positions!! 

In the resulting Lorentzian model, the normal state 
E is purely an imaginary constant, and a 2 F is a mode 
at zero energy. In the superconducting state, this mode 
gets pushed back to some energy within 2 A. This model 
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is artificial in the sense that all the self-energy is being 
generated by the mode. That is why we went through the 
above discussion motivating the mode more properly as 
a rearrangement of a 2 F due to the superconducting gap. 
In practice, though, the results are very similar to the 
MFL plus mode model, and has the further advantage 
of having the several parameters of that model collapse 
to just the mode strength (Ti) and mode position (Hq) 
of the Lorentzian model. Moreover, analytic results can 
still be obtained for E when the superconducting density 
of states for the k — q line of Fig. 1 is taken into account. 
The result is 

- JmE(w) = IW(|w|) + TiN(\u\ - ft ), \w\ > £! + A 

= iyv(|w|), a< m <n + A 

= 0, |w|<A (2) 



where N(u>) = lj/^/oj 2 — A 2 is the BCS density of states, 
and 



7ri?eE(w) = r N(-uj) In \-u+ ^Juj 2 - A 2 |/A 



+TiN(ft - w) In |fi - w + - Q ) 2 - A 2 |/A 

-{cj^-u;} (3) 

where it has again been assumed that A is a real con- 
stant in frequency. An s-wave density of states has been 
used to obtain an analytic result. A d-wave density of 
states will not be that different. The advantage of an 
analytic result is that it is useful when having to take 
spectra and convolve with resolution to compare to ex- 
periment. Our results are not very sensitive to Tq (30 
meV), included again to damp the quasiparticle peak. 
(A more realistic damping of the peak would require 
making A complex.) We use the same set of parameters 
for all k (Ti=200 meV), and therefore assume a d-wave 
gap Ak = A max (cos(k x a) — cos(fc y a))/2 in Eqs. 1-3 with 
A max = 32 meV. The best agreement with experiment is 
found by choosing the mode energy Slo = l-3Ak, so that 
the spectral dip for (n, 0) is at 2.3A max . 

The resulting real and imaginary parts of £ at (tt, 0) 
are shown in Fig. 3a. Note the singular behaviors at 
A due to the Tq term and at £Iq + A due to the Ti 
term. In both cases, step drops in JmS would also give 
singularities in ReY>. The advantage of peaks in JmS 
(due to the density of states) is that it makes the dip 
deeper in better agreement with experiment. In Fig. 3b 
and 3c, we show a comparison of the resulting spectral 
function (convolved with the experimental energy and 
momentum resolution) to experimental data at (tt, 0) for 
both wide and narrow energy scans, where a step edge 
background with a gap of A is added to the calculated 
spectrum as discussed above. The resulting agreement is 
excellent. 
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FIG. 3. (a) JmS and ReZ at (tt,0) from Eqs. 2-3 
(ri=200meV, r =30meV, A=32meV, n =1.3A). Compari- 
son of the data at (tt, 0) for (b) wide and (c) narrow energy 
scans with calculations based on Eqs. 1-3, with an added step 
edge background contribution. 

To better appreciate these results, the positions of the 
sharp peak and the higher binding energy hump obtained 
from the calculations are plotted relative to the normal 
state binding energy et along the (0, 0) — (tt, 0) direction, 
and compared to those obtained from the experimental 
data of Ref. |] in Fig. 4. This plot is very similar to 
that obtained for electrons interacting with an Einstein 
modep. The calculations reproduce the dispersionless na- 
ture of the low frequency peak, as well as its lack of visi- 
bility for k vectors close enough to (0,0). The dispersion- 
less behavior is due to several factors: (1) the weakness of 
the dispersion ek near (tt, 0), (2) the lowering of f2 ~ Ak 
as one moves towards (0, 0), and (3) the influence of both 
ReYi and A. The last is a new effect worth commenting 
on. The real part of the self-energy implies a mass en- 
hancement (Z > 1) in the superconducting state relative 
to the normal state, which acts to push spectral weight 
towards Ep . On the other hand, A itself pushes spectral 
weight away from Ep. Thus the dispersion is dramati- 
cally flattened relative to the normal state. 
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FIG. 4. Positions (meV) of the sharp peak and the broad 
hump in the superconducting state versus normal state peak 
position. Solid points connected by a dashed line are the 
experimental data, the solid lines are obtained from the cal- 
culations, and the dotted line represents the normal state dis- 
persion. 

In the above calculations, it was assumed that the 
mode frequency was proportional to Ak- This was the 
easiest way we found to properly simulate the loss of 
the experimental low frequency peak as one disperses to- 
wards (0, 0). In reality, Flo is a function of q in the dia- 
gram of Fig. 1, not k. Moreover, it was our assumption 
of independence of f2n on q that allowed us to obtain a 
step drop in JmS, leading to the spectral dip. Without 
some microscopic theory, only qualitative observations 
can be made at this stage concerning the true depen- 
dence of E on momentum kE£l. Assuming an artificial 
limit where only q = Q = (w, ir) contributes, we would 
replace TiN(uj + fi ) in Eq. 2 by g\_ k+Q A k+ Q(uj + fi ) 
(for lu < 0) where g is the interaction vertex. Using 
a quasiparticle pole approximation for A when solving 
Eqs. 2-3, this would imply a dip in the spectrum at 
\oj\ = E k+Q + tt where £f = e 2 k /(ReZ k ) 2 + A 2 k , and a 
persistent low frequency peak if Z is large enough. The 
coupling of k and k + Q in the self-energy equations also 
implies that if a low frequency peak exists for k, then 
one also exists for k + Q. This_is just the effect observed 
in the data along (ir, 0)-(tt, ir)U, in that a low frequency 
peak exists for about the same momentum range as that 
along (ir, 0)-(0, 0). It remains to be seen whether such 
simple momentum dependent models give as good a fit 
to the spectra as the dispersionless model presented here. 

Finally, it is interesting to note that the mode en- 
ergy we infer from the data is 41meV, equivalent (prob- 
ably fortuitously) to a resonant made energy observed in 
YBCO by neutron scattering datao at Q = (7r, it). The 
models proposed for this mode are similar to the model 
discussed in this papenl3. So far, neutron scattering data 
on Bi2212 have yet to see a similar structuretLl3, although 
these experiments were done on a rod of aligned small 
crystals. Our results here would imply that such experi- 



ments on large single crystals would be of interest. 

In conclusion, we have shown that a simple model of 
an electron interacting with a collective mode in the su- 
perconducting state gives a quantitative description of 
the unusual spectral lineshape seen by ARPES data in 
the superconducting state of Bi2212. This implies that 
electron-electron scattering plays a dominant role in high 
temperature superconductors, and is in support of an 
electron-electron origin for the pairing. 
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